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INTRODUCTION 

/X« 


This  report  is_-the  seventh- bimonthiv  progress-Teport ' fflverlng 
,  A 

.SRI-  ImrcntaTiunai-^  current  study  of  the  response  of  adaption  kits  in, 
accidental  side  impacts.  The  primary  goal  of  the  project  is  the 
further  development  of  finite  element  models  for  a  generic  missile 
structure  so  that  adaption  kit  response,  particularly  component 
accelerations,  can  be  predicted  for  side  impacts  of  the  missile.  An 
additional  goal  is  experimental  determination  of  the  impact  response 
of  a  missile  structure  that  contains  a  hard  link  safe-arm  device. 


The  work  planned  far-frhia  -eeotraet  consists  of  four  tasks.  In  ^ 
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Task  1,  we  are  furehaj  developing  finite  element  models  for  substructures 
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and  checking  the  models'  predictions  against  experimental  data. 

Task  2,  w^re^  the  models^to  ^iVulVtt  the  response  of  complete 
structures,  again  by  comparing  the  predictions  with  experimental  results 
insofar  as  feasible. 


In  Task  3,  -we  are  transferring  the  analytical 
-  t  “e  ,  V..-AX 


capabilities  developed  in  Tasks  1  and  2xby  implementing  the  SUPER  code 
at  ARRADCOM.  In  Task  4,  wr  are-pe*£ermtwg  impact  experiments^on  scale 
models  of  a  specific  prototype  missile  structure  to  help  plan  some  full- 
scale  tests,  fafor  a  more  detailed  description  of  Tasks  1  through  3,  refer 
to  Bimonthly  Progress  Report  No.  1,  July  1978.  Task  4  is  described  fully 
in  Bimonthly  Progress  Report  No.  3,  November  1978. 

During  this  reporting  period,  progress  was  made  in  Tasks  1,  2,  and  3. 
The  estimated  percentage  of  the  work  completed  is  shown  in  Figure  1.  The 
remainder  of  the  report  describes  the  technical  progress  and  future  plans 
for  each  task  and  indicates  the  financial  status  of  the  contract. 
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PROGRESS 

Task  1;  Development  of  Finite  Element  Models  for  Structures 

Improved  models  for  the  ring,  bolts,  and  plate  are  being  developed. 
The  ring  model  has  been  improved  by  including  a  better  representation  of 
the  shell  inertia.  The  bolt  model  will  include  unloading  and  reloading 
The  plate  model  has  been  extended  to  include  bending. 

Ring  Model 

Four  comparisons  have  been  made  between  the  calculated  and 
measured  response  of  rings  with  rectangular  cross  sections  (flat  ring) 
or  with  L-shaped  cross  sections  (L-ring);  these  comparisons  will  be  pre¬ 
sented  fully  in  the  next  report.  Comparisons  were  made  for  four  cases: 

(1)  a  flat  ring  at  22  ft/sec 

(2)  an  L-ring  at  17  ft/sec 

(3)  an  L-ring  at  34  ft/sec 

(4)  an  L-ring  with  the  shell  at  17  ft/sec. 

In  each  comparison  the  agreement  between  the  analysis  and  the  experi 
ment  is  good.  In  the  investigation  of  the  ring  model,  it  was  discovered 
that  the  very  high  frequency  response  predicted  previously  is  not  excited 
in  the  model  when  shell  effects  are  included.  Comparison  between  the  pre 
dieted  and  measured  ring  responses  will  be  included  in  the  next  report. 

Bolt  Model 

It  is  anticipated  that  the  current  bolt  model1  will  be  inappro¬ 
priate  for  three-dimensional  calculation  due  to  the  absence  of  an  axial 
force.  A  new  bolt  model  is  being  formulated  for  use  in  plate  bending 
calculations.  The  new  bolt  model  will  be  described  in  a  future  report. 
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Plate  Model 


The  triangular  plate/shell  element  used  in  the  SUPER  code  is  based 
on  the  work  of  ‘larchertas  and  Belytschko-  for  the  explicit  formulation 
and  Kulak  and  Belytschko3  for  the  implicit  formulation.  This  formulation 
uses  a  co-rotational  or  convected  coordinate  system.  These  coordinates 
facilitate  the  formulation  of  element-related  quantities;  e.g.,  they 
allow  simple  strain-nodal  displacement  as  well  as  nodal  force-stress 
relations.  Argyris  et  al.u  used  co-rotational  coordinates  in  static 
problems . 

The  convected  coordinate  formulation  uses  three  coordinate  systems: 
(1)  a  fixed  global  coordinate  system  designated  by  (x,  y,  2);  (2)  a  body 
coordinate  system  (x,  y,  z)  for  each  node  that  rotates  with  the  node  and 
coincides  with  the  principal  axes  of  the  mass  moment  of  inertia  of  the 
node;  (3)  a  convected  coordinate  system  (x,  y,  z)  originating  at  node  1 
of  the  element,  and  following  the  element  so  that  the  (x,  y)  plane  always 
passes  through  the  three  nodes  of  the  element  and  the  x-axis  bisects  the 
current  angle  at  node  1  (see  Figure  2).  Thus,  the  (x,  9>  z)  coordinate 
system  approximates  the  rigid  body  motion  of  the  element.  The  plate/shell 
element  can  be  subjected  to  linear  in-plane  and  cubic  transverse  dis¬ 
placements.  The  displacement  of  each  element  is  decomposed  into  rigid 
body  motion  and  deformation.  The  rigid  body  motion  is  defined  by  the 
motion  of  the  convected  coordinate  system  (x,  y,  z)  .  The  deformation 
displacements  are  the  displacements  relative  to  the  convected  coordinates. 
Thus,  the  deformation  displacements  may  be  obtained  by  subtracting  the 
rigid  body  displacements  from  the  total  displacements  or  by  characterizing 
them  by  modes  that  are  independent  of  the  rigid  body  motion. 

The  in-plane  deformation  is  calculated  from  the  elongation  of  the 
element  sides,  and  the  transverse  deformation  is  related  to  nodal  ro¬ 
tations  relative  to  the  (x,  y)  plane.  These  nodal  rotations  are  obtained 
by  the  vector  cross  product  of  the  current  element  unit  normal  and  the 
original  element  unit  normal  rotated  through  the  element's  rigid  body 
motion. 


Element  strains  are  calculated  by  differentiating  the  deformation 
displacements  and  by  using  the  Kirchoff  assumption  that  normals  to  the 
midplane  remain  straight  and  normal.  Currently,  only  the  Von  Hises 
yield  criterion  with  isotropic  hardening  can  be  used  to  calculate  element 
stresses,  although  other  constitutive  relations  can  be  incorporated. 

As  in  the  case  with  mose  plate/shell  bending  elements,  the  element 
formulation  described  by  Belytschko  and  co-workers2’3  is  a  nonconforming 
element.  The  nonconformity  for  this  element  is  that  it  violates  the 
continuity  of  slope  conditions  and  therefore  only  approximates  the  prin¬ 
ciple  of  minimization  of  potential  energy.  The  continuity  of  slope 
3w/3n  can  be  satisfied  at  the  nodes  but  not  along  the  entire  edge  of 
the  element  because  for  a  cubic  displacement  field,  the  slope  3w/3n  is 
is  quadratic  along  an  edge.  The  following  example  illustrates  this 
point.  Consider  the  triangular  bending  element  shown  in  Figure  3  with 
an  element  coordinate  system  along  side  1-2  denoted  by  x  and  y.  Along 
the  side  1-2  we  have 


3w  3w  3w  3x  3w  3y 
3n  *  3y  3x  3y  3y  3y 


-sin9  +  cosQ 

3x  3y 


or  w  .  ■  a  w  where  a  is  the  usual  tensor  transformation.  For  a 
» 1  lj  » J  ij 

cubic  displacement  function 

w  -  aL  +  a2x  +  a3y  +  a4x2  +  a5xy  +  a&y2  +  a?x3  +  agx2y  +  agxy2  +  a1Qy3 

This  function  would  have  the  same  form  in  barred  (-)  coordinates. 

Thus, 


3w 

3y 


a3  +  a5  x  +  2afiy  +  agx“2  +  2agyx  +  3aI0y“- 


Along  the  edge  1-2,  y  *  0;  hence 


3w 

3y 


I3  +  a5  x  +  I8x2 
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and  we  see  that  the  slope  3w/3n  ■  3w/3y  is  a  quadratic  function  along 
the  edge.  Since  we  have  only  two  conditions  along  the  edge,  i.e., 

3w/3n  at  node  1  and  2,  we  cannot  determine  all  of  the  constants  uniquely. 
The  most  common  remedy  for  this  problem  is  to  omit  the  xv  term  (i.e.,  set 
a^«0)  in  the  cubic  equation  for  displacement.  Consequently  the  element 
is  nonconforming  since  its  displacement  field  is  represented  by  an  in¬ 
complete  cubic  polynomial. 

Although  the  bending  element  formation  is  nonconforming,  this  does 
not  mean  that  solutions  using  these  elements  will  diverge  from  the 
correct  response.  It  has  been  shown  by  Irons  and  co-workers-  that 
the  simple  triangular  element  converges  to  an  exact  solution  if  the  mesh 
is  generated  by  sets  of  parallel  and  equally  spaced  lines.  Also,  for 
practical  engineering  purposes,  in  most  cases  the  accuracy  obtained  by 
the  nonconforming  triangle  is  adequate.  Indeed,  at  most  practical  sub¬ 
divisions,  it  gives  results  superior  to  those  attainable  with  equivalent 
conforming  triangles  (see  Reference  5). 

To  determine  the  accuracy  of  the  plate/shell  bending  element  for 
the  adaption  kit  response  problem,  we  solved  two  plate-bending  problems: 
static  and  dynamic  loading  of  a  simply  supported  circular  plate  under 
a  uniformly  distributed  load.  The  mesh  configuration  used  for  these 
calculations  is  shown  in  Figure  4.  This  nonuniform  discretization  was 
chosen  because  it  is  typical  of  the  mesh  configurations  used  in  the 
adaption  kit  response  simulations.  Figure  5  compares  the  static  response 
predicted  by  the  SUPER  code  with  the  analytical  solution  obtained  from 
Reference  6.  This  comparison  indicates  that  the  code  over-predicts  the 
displacements  by  about  4%.  Similarly,  Figure  6  shows  the  results  of 
a  dynamic  calculation  and  the  analytical  solution  obtained  from 
Reference  7.  In  this  case,  the  code  is  shown  to  underestimate  the  de¬ 
flections  by  an  average  of  8%.  Overall,  these  results  agree  with  those 
reported  by  other  users  of  the  nonconforming  elements;  for  example, 

R.  Szilard  reported  differences  of  6%  between  computed  and  analytical 
results.  Of  course,  a  more  meaningful  comparison  will  be  made  between 
the  computed  and  measured  responses  of  the  generic  missile  structure. 
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Task  2:  Calculation  of  the  Response  of  Complete  Structure 

An  experiment  with  a  complete  structure  at  an  impact  velocity  of 
80  ft/sec  was  conducted  to  measure  bending  response  and  to  determine 
whether  any  new  response  mechanisms  occur  at  high  speed.  Six  strain 
gages  were  mounted  on  the  plate  to  measure  top  and  bottom  surface  strain 
in  front  of  the  AK  mass,  behind  the  AK  mass,  and  to  the  side  of  the  AK 
mass.  Two  accelerometers  were  mounted  on  the  AK  mass  to  measure  radial 
and  axial  acceleration. 

The  results  of  the  experiment  are  plotted  in  Figures  7  through  10. 

The  measurements  show  that  large  bending  strains  occur  during  the  early 
part  of  the  response,  the  radial  acceleration  is  limited  by  the  mounting 
bolts,  and  the  axial  acceleration  is  significant.  Further  analysis  of 
this  experiment  is  in  progress. 

Task  3:  Implementation  of  SUPER  at  ARRADCOM 

A  tape  of  the  SUPER  code  with  example  data  has  been  sent  to  ARRADCOM. 
A  tape  of  the  final  version  of  the  code  will  be  sent  at  the  completion  of 
the  contract. 

Task  4:  Model  Tests  of  Missile  with  Hard  Link 
This  task  is  complete. 
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FINANCIAL  STATUS 

As  of  2  July  1979,  $63,457  has  been  spent  on  labor  (1282  supervisory 
and  professional  hours,  729  technical  and  clerical  hours),  and  $8,216 
has  been  spent  on  materials  and  services.  Of  the  total  contract  funds 
$93,441,  the  balance  remaining  is  $21,768. 
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Reporting  Period:  3  May  -  2  July  1979 


Hours 


Total  hours  expended  to  date: 


Supervisory  and  Sr.  Professional  Personnel  267 
Professional  Personnel  1015 
Technician  729 

Cumulative  total  hours  to  date  2011 
Percent  of  total  hours  expended  to  date  83% 


Funds 

Funds  expended  during  the  reporting  period  $14,245 

Funds  expended  to  date  71,673 

Percent  of  total  funds  expended  to  date  77% 


Work 

Percent  of  work  completed  during  the  reporting 

period  10% 

Percent  of  work  completed  to  date  $0% 
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{ 


FIGURE 


MESH  FOR  SIMPLY  SUPPORTED  PLATE 
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FIGURE  5  STATIC  PLATE  BENDING  COMPARISON 
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